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It  is  shown  that  cross  field  temperature  gradients  can  drive  ion 
acoustic  waves  unstable  in  a laser  produced  plasma.  The  principle 
effect  of  the  instability  is  an  enhancement  of  the  cross  field  thermal 
conductivity.  It  is  shown  how  this  effect  can  be  modeled  in  laser 
fusion  hydrocodes. 
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CROSS  FIELD  THERMAL  TRANSPORT 
DUE  TO  ION  ACOUSTIC  WAVES  IN  MAGNETIZED 
LASER  PLASMAS 

I.  INTRODUCTION 

IC  Is  now  well  established  that  magnetic  fields  can  be  spontane- 
ously generated  In  laser  produced  plasmas.  Recent  measurements1  by 
Faraday  rotation  have  Indicated  fields  larger  than  one  megagauss  In 
the  underdense  plasma.  There  are  now  many  mechanisms  by  which  ma>  a- 
tic  fields  can  be  produced.  These  include  the  7n  x TT  term  In  the 
equation  for  B,2*3  the  thermo-electric  term,2"4  the  divergence  of 
momentum  flux  of  laser  light,5*®  currents  produced  by  energetic  elec- 
trons,' and  Impurity  seeding  and  even  currents  generated  by  turbulence 
in  an  unmagnetized  plasma.9 

Although  the  magnetic  field  Is  large,  the  plasmas  are  so  dense, 
so  hot  and  are  flowing  so  rapidly  chat  the  magnetic  fields  observed 
and  predicted  are  not  nearly  large  enough  to  have  any  effect  on  the 
dynamics.  The  principle  way  In  which  the  effect  of  the  magnetic  field 
makes  Itself  felt  Is  by  greatly  reducing  the  cross  field  electron 
thermal  conduction.  In  order  for  the  thermal  conduction  to  be  reduced, 
one  need  only  have  u>  t » 1,  where  is  the  electron  cyclotron 
frequency  and  t is  the  electron  collision  time.  This  condition  is 
very  easy  to  satisfy  even  for  magnetic  fields  too  small  to  affect  the 
dynamics. 
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Thus  one  expects  Che  cross  field  thermal  energy  flux  In  a magnet- 
ized plasma  to  be  much  less  than  Che  flux  In  an  unmagnetlzed  plasma. 
This  Inhibition  of  energy  flux  would  In  turn  lead  to  much  steeper 
electron  temperature  gradients  than  what  would  be  expected  for  an  un- 
magnetlzed  plasma.  One  could  then  ask  whether  the  steep  temperature 
gradients  so  produced  could  lead  to  some  sort  of  Instability  which 
would  enhance  the  electron  thermal  flux  and  reduce  the  temperature 
gradient.  Experience  with  magnetically  confined  plasmas  certainly 
indicates  chat  this  will  be  the  case.  For  instance  In  tokamaks, 
electron  temperature  gradients  drive  trapped  particle  Instabilities*' 
which  are  thought  to  greatly  enhance  the  electron  thermal  conduction. 
In  any  case  It  Is  a fact  chat  Che  electron  thermal  conduction  In 
tokamaks  is  enhanced  over  Its  classical  value  by  many  orders  of 
magnitude.11 

In  this  paper,  we  show  that  ion  acoustic  waves  which  propagate 
perpendicular  to  both  the  magnetic  field  and  the  gradients  can  In 
fact  be  driven  unstable.  Section  II  derives  Che  equilibrium.  Section 
III  calculates  the  linear  theory  for  Ion  acoustic  waves  In  this 
equilibrium.  Section  IV  calculates  the  electron  energy  flux  In  a 
given  spectrum  of  Ion  acoustic  fluctuations.  It  Is  shown  there  that 
the  electron  thermal  conduction  scales  as  |e©/T|2.  We  also 
calculate  how  small  |eo/T|2  must  be  In  order  that  the  magnetic  field, 
rather  than  the  turbulence  in  an  unmagnetlzed  plasma,  Is  primarily 
responsible  for  Inhibiting  the  electron  thermal  flux.  Finally  in 
Section  V,  we  show  how  to  Include  the  effects  of  enhanced  cross  field 
transport  In  a fluid  code. 
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II.  THE  EQUILIBRIUM 

The  basic  configuration  of  a one  dimensional  laser  produced  plasma 

is  shown  In  Fig.  1.  The  laser  Is  to  the  left  at  x =-*.  The  flow 

velocity  Is  In  the  negative  x direction,  the  temperature  gradient  Is  In 

the  negative  x direction  and  the  density  gradient  Is  In  the  positive  x 

direction.  Also,  we  assume  the  existence  of  a magnetic  field  B In  the 

positive  Z direction  (out  of  the  plane  of  the  paper).  The  magnetic 

field  is  assumed  to  be  strong  enough  to  affect  the  electron  transport, 

i.e.,  uu  t » 1 where  w Is  the  electron  cyclotron  frequency  and 
ce  e ce 

t Is  the  electron  collision  time.  However,  It  Is  also  assumed  to  be 
e 

sufficiently  weak,  that  It  has  no  effect  on  the  dynamics,  l.e., 

nMu^  » B2/Vr  » ou  and  P » l where  n,  M,  u and  u>  are  reapectlvely 

pe  ce  pe 

the  number  density,  Ion  mass,  flow  velocity  and  electron  plasma  frequency 

and  6 ■ LrrnT/B2.  Also,  the  Ions  are  assusted  to  be  unaagnetlzed,  chat  Is 

0 ^ » L where  la  the  Ion  larmor  radius  and  L Is  a macroscopic  scale 

length.  On  the  other  hand  the  electrons  are  assumed  to  be  strongly 

magnetized,  or  « L.  Therefore,  in  order  for  the  electrons  to  flow 

In  the  negative  x direction  there  must  be  an  elactrlc  field  In  the 

negative  y direction  with  magnitude  B u/c.  The  above  conditions  are 

generally  well  satisfied  in  a laser  produced  plasma. 

In  order  to  calculate  whether  such  a plasma  is  stable,  it  is 

essential  to  consider  the  effect  of  collisions  on  the  equilibrium. 

The  equilibrium  is  governed  by  a flow  velocity  to  the  left,  and  an 

electron  thermal  conduction  to  the  right.  This  thermal  energy  flux 
dX 

is  equal  to  Q = -K  — - where  K is  the  thermal  conductivity.  Classical 
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kinetic  theory  gives  the  result:1' 

K = k.66  n T /our  r (l) 

e c e 

If  the  laser  energy  flux  is  I and  the  fractional  absorbtion  is  Or,  then 
conservation  of  energy  flux  (assuming  for  simplicity  that  there  is  no 
change  in  fluid  energy  flux)  across  the  region  of  absorption  gives 

Q = al.  (2) 

Thus  thermal  energy  flux,  and  therefore  a non-zero  collision  frequency, 
is  inherent  in  the  very  nature  of  the  equilibrium  of  a laser  produced 
plasma.  The  Appendix  discusses  further  the  lack  of  Vlasov  equilibrium 
for  the  configuration  shown  in  Fig.  1. 

As  we  will  see,  cross  field  thermal  conduction  in  the  x direction 
is  always  accompanied  by  an  additional  thermal  flux  in  the  y direction 
(along  iso  therms).  The  magnitude  of  this  thermal  flux  in  the  y direc- 
tion is  independent  of  collisions  (for  u>ce  « 1)  at  given  tempera- 
ture gradient.  This  thermal  flux  perpendicular  to  both  field  and 
gradients  is  well  known  in  classical  kinetic  theory.13  It  is  analogous 
to  diamagnetic  currents  in  an  Isothermal  plasma,  which  are  perpendic- 
ular to  both  the  magnetic  field  and  density  gradient. 

In  the  next  section,  it  is  shown  that  this  thermal  flux  in  the  y 
direction  gives  rise  to  the  possibility  of  ion  acoustic  instability, 
with  wave  vector  in  the  y direction.  Since  this  energy  flux  in  the 
y direction  is  nearly  independent  of  collision  frequency,  we  will  use 
the  simplest  possible  collision  model,  a Krook  collision  term.  The 
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steady  state  Vlasov  equation  becomes: 


V - - E £ E r—  - ~ - X B * |4  ■ (3 

x dx  m x dV  m y oV  m c w oV 
x ' y » 

where  f°  Is  a Maxwellian  distribution  with  the  local  density  n,  flow 
velocity  u and  temperature  T , which  are  x dependent.  If  we  make  the 
simplifying  assumption  that  v Is  Independent  of  V,  the  collision  term 
conserves  particles,  momentum  and  enargy.  The  quantity  E^  cannot  be 
specified  Independently.  It  is  determined  by  the  x component  of  the 
electron  momentum  equation. 


0 ' --  EX  - 51  " V 


where  we  have  assumed  u2  « T^/m.  Then  assuming  that  the  first  two 
terms  on  the  left  hand  side  of  Eq.  (3)  are  small,  one  can  solve  for 
the  perturbed  distribution  function  by  standard  means.14  Assuming 


fO  n(x) 

(2tt  T (x)/m)a 


exp  - 


m (V  - ui  Y 
2 Te(x) 


f = f + 6 f , 


Eq.  (3)  can  then  be  manipulated,14  by  using  Eq.  (k) , into  an  equation 


for  6f 


“c  d6  5f  + v6f  = -(vx-u> 


. dT  f m (V  - ui  )2  \ 


by  making  use  of  the  transformation 


V - u = V 0r9 
x i 


(8) 

V = V s in  9. 

y j- 


We  have  retained  on  the  right  hand  side  of  Eq.  (7)  only  terms  which 
give  rise  to  a nonzero  energy  flux.  Solving  Eq.  (7)  we  find 


f * fQ+«f  = <1 


T , -<S-Uix>2  *1  1 « /*<Vx-«)»  "cVvU 

'I  Te  ' 2J  Te  dx  \ V4  ^ yj 


fo  <9) 


The  first  term  in  the  parentheses  of  Eq.  (9)  gives  rise  to  an  energy 
flux  down  the  temperature  gradient.  If  this  energy  flux  is  imposed  as 
a boundary  condition,  then  the  temperature  gradient  is  thereby  spec- 
ified. The  second  term  in  the  parentheses  of  Eq.  (9)  gives  rise  to  an 
energy  flux  in  the  positive  y direction.  In  the  limit  of  v « a>c, 
this  energy  flux  in  the  y direction  is  Independent  of  colllsionallty 
and  it  depends  only  on  the  temperature  gradient.10  Notice  that  the 
energy  flux  in  the  positive  y direction  is  carried  by  a flux  of 
energetic  electrons  moving  in  the  positive  y direction  and  that  this 
flux  is  balanced  by  a return  current  of  slower  moving  electrons  in  the 
me gative  y direction.  In  the  next  section,  we  will  show  that  this 
return  current  in  the  negative  y direction  can  drive  instable  ion 
acoustic  waves  which  propagate  in  the  negative  y direction. 

Ill  LINEAR  THEORY 

The  calculation  of  the  perturbed  electron  distribution  function, 
given  the  unperturbed  distribution  function  in  Eq.  (9)  can  be  done 
with  standard  techniques. 15  The  result  is 
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Acv  \r  . kc  dT«4  m(.v  - UixV  s\l 

— JJ 


-lx  -<*> 

ce 


1°  (10) 


assuming  ChaC  k is  in  Che  y direction  Che  fluctuating  potential  is, 


o(y,t)  = ® exp  i (ky  - uut)  c.c. 


and  cu  » v. 

ce 

At  this  point,  we  make  use  of  the  fact  that  w » u»  and  kV  » w 

ce  e c 

where  = (T^/m)  . A great  deal  of  analytic  work  and  also  numerical 
simulation  has  shown  that  In  this  limit  the  magnetic  field  has  no  effect 
on  the  microscopic  dynamics  of  the  Instability.  However,  the  dT/dx 
term  in  Eq.  ( 10 ) is  unchanged.  (This  point  will  be  discussed  more 
fully  in  the  next  section.) 

Crudely,  one  can  envision  the  situation  as  follows.  As  expressed 

in  Eq.  (10),  is  a summation  over  a large  number  of  resonances.  It 

has  been  shown  that  any  physical  mechanism  which  broadens  the  cyclotron 

resonances  by  the  cyclotron  frequency15  has  the  effect  of  converting 

the  magnetized  dispersion  relation  into  an  unmagnetized  dispersion 

relation.  The  broadening  mechanism  may  either  be  collisions,  resonence 

broadening1^  or  cyclotron  trapping.17  Furthermore,  if  k^  is  not  exactly 

zero  but  has  some  small  value  (k  V > u>  ) it  is  now  well  established 

z e ce 

that  even  for  linear  theory  in  a collisionless  plasma,  the  cyclotron 
resonances  are  washed  out.13 

For  the  case  in  which  resonance  boradening  governs  the  resonance 
width,  it  has  been  shown  that  if  k \D#  ~ 0.5,  (the  case  we  will  examine 
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here>  Chen  Che  Individual  cyclocron  resonances  are  washed  ouC  if 


*r  19  /*  ) ',s. 


ce  pe 


Equacion  (12)  above  shows  chat  very  small  values  of  ec/T  are  sufficienc 
Co  eliminace  Che  individual  cyclocron  resonances.  Specifically,  Che 
values  of  eo/T  needed  Co  wash  ouC  individual  resonances  are  much  less 
Chan  chac  needed  Co  wipe  ouc  Che  effecC  of  Che  magnetic  field  on  Che 
equilibrium  recura  currenc  in  che  y dlrecclon.  In  che  nexc  secCion,  lc 
is  shown  chac  che  magnecic  field  dominaces  Che  equilibrium  currenc  and 
cransporc  as  long  as 


Thus,  chere  is  a large  range  of  — for  which  che  microscopic 

e 

picture  is  governed  by  waves  in  an  unmagnetized  plasma,  buC  for  which 
Che  macroscopic  currents  and  CransporC  are  dominated  by  Che  magnetic 
field.  This  Chen  is  similar  Co  che  case  of  a resistive  shock  in  a 
dense  plasma  >:>  where  the  magnetic  field  governs  the  currents, 

but  where  the  instability  responsible  for  anomalous  transport  is  an  ion 
acoustic  instability  of  an  unmagnetized  plasma. le~20 

Than  averaging  over  cyclotron  resonances,  the  expression  for 


f becomes 
a 


f = nf°  2* 
e T 

e 


r + 

L1+" 


. dT  / m 

— —*  ( i - 

eB  dx  y 


g - 0 


k*V-<D 


Assuming  Chat  che  ions  are  cold  and  unmagnetized , and  k2\2  « 1,  the 
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dispersion  relation  for  waves  is  obtained  by  n^  = n^,  or 


(IT 

-Ei  = —££_ 

x-  k-V* 


- AT  / m(V  " \ 

f .*  H £(»  — T r--  .?) 

1*  f!  ,B  e “''•I 

J y J 


(15) 


dT 


C g 

As  long  as  — « V^,  the  Integral  term  In  the  square  brackets 

is  much  smaller  than  unity,  so  to  lowest  order  the  waves  are  just  Ion 

acoustic  waves  which  propagate  with  phase  speed  V = (T  /M)^.  Making 

8 0 

a resonant  approximation  to  (kV^  - tu)”1,  we  find  that  the  wave  fre- 
quency is  given  roughly  by 


x * I k|v 


* 1 \*  * J 

l • i("/e>  — tct 


(■  - stS) 


(l  6) 


dT 


In  the  configuration  shown  in  Fig.  1,  ^ < 0,  so  that  if  k < 0, 
the  ion  acoustic  wave  is  unstable.  This  corresponds  to  a wave  propa- 
gating in  the  negative  y direction,  or  parallel  to  che  return  current 
in  the  y direction.  For  the  case  of  cold  ions  which  we  have  considered 
here,  the  condition  for  instability  is  given  by 


5/2  dr  -k > 


(17) 


This  condition  is  easily  satisfied  in  a laser  produced  plasma.  For 
Instance  if  T = 5 Rev  and  B * 1 Magagauss,  a hydrogen  plasma  is  un- 
stable if  the  temperature  gradient  scale  length  is  less  than  about 
75  U.  If  the  ion  temperature  is  not  zero,  one  can  subtract  from  the 
growth  rate  the  ion  Landau  damping  decrement. 
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IV  ANOMALOUS  THERMAL  CONDUCTION 


In  this  section,  we  will  calculate  the  anomalous  thermal  conduction 
arising  from  the  ion  acoustic  instabilities  driven  by  cross  field  tem- 
perature gradients.  As  discussed  in  Section  II  and  III,  the  negative 
temperature  gradient  in  the  x direction  gives  rise  to  a current  of  low 
velocity  particles  in  the  negative  y direction.  This  drives  an  ion 
acoustic  Instability  which  also  propagates  in  the  negative  y direction. 

One  consequence  of  the  instability  is  that  a force  is  exerted  on 
the  electrons  opposite  to  the  direction  of  the  drift,  or  in  the  positive 
y direction.  This  force  acts  on  the  electrons  like  an  equivalent  elec* 
trie  field  in  the  negative  y direction.  This  equivalent  electric  field 
crossed  into  the  magnetic  field  in  the  positive  z direction  gives  rise 
to  a net  electron  drift  in  the  negative  x direction  i.e. , up  the 
temperature  gradient. 

This  hydrodynamic  motion  of  the  electrons  of  course  is  not  real. 

In  a one  dimensional  laser  produced  plasma  characterized  by  <u  » uu^, 

S » 1,  and  nMu2»  the  magnetic  field  has  no  effect  on  the  hydro- 

dynamic  motion.  This  is  controlled  only  by  the  ion  motion  and  the 
quasi-neutrality  condition.  (In  a two  or  three  dimensional  plasma, 
the  magnetic  field  can  play  a role  in  the  electron  motion  as  long  as 
7*ue  ~ 7*V>  Therefore  the  electric  field  in  the  y direction,  which 
controls  the  x motion  of  the  electrons  must  also  be  modified  by  the 
instability. 

In  the  reference  frame  moving  with  speed  cE^/B,  one  can  see  how 
the  instability  gives  rise  to  an  anomalous  thermal  flux.  In  Fig.  2a, 


I 
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the  dotted  lines  are  the  contours  of  constant  f°(V)  in  the  (V  -u) , V 

e x y 

plane  in  the  absence  of  instability.  The  instability  exerts  an  average 
force  on  the  electrons  in  the  positive  y direction  which  causes  an 
average  electron  drift  in  the  negative  x direction.  However  this  force 
is  a decreasing  function  of  the  magnitude  of  the  particle  velocity.'  ‘ 
Thus  the  lower  velocity  particles  have  a greater  drift  velocity  in  the 
negative  x direction  than  do  the  higher  velocity  particles.  Hence  in 
the  presence  of  Instability,  the  contours  of  constant  f(V)  are  the 
solid  lines  in  Fig.  2a. 

Transforming  to  the  reference  frame  in  which  there  is  no  net  drift 

in  the  x direction,  the  contours  of  constant  f(V  ) are  as  shown  in 

€ 

Fig.  2b.  Notice  that  the  lower  velocity  particles  have  a drift  to  the 
left  (up  the  temperature  gradient),  whereas  the  higher  velocity  parti- 
cles have  a drift  to  the  right  (down  the  temperature  gradient).  Since 
(V  ) = 0,  this  distorted  distribution  function  has  an  energy  flux  down 
the  temperature  gradient. 

As  long  as  u2  « T^/m,  a condition  easily  satisfied  in  a laser 
produced  plasma,  the  electron  energy  flux  in  Fig.  2b  is  just  equal  to 
the  energy  flux  in  Fig.  2a  minus  5/2{nuVTe,  where^nu^is  the  instability 
induced  particle  flux.  The  quantity  is  of  course  just  that 

part  of  the  total  energy  flux  which  is  convected  with  the  fluid  motion 
(for  u2  « T^/m). 

We  will  now  calculate^uj^  and  W,  the  particle  and  energy  flux 
associated  with  the  instability.  If  we  make  transformation  given  in 
Eq.  (6),  the  quasl-linear  equation  for  f(V)  in  a time  independent 
plasma  is 

1 1 


2 


(18) 


+ w 


'll 

C *0 


k 

«< 


«c(k) 

m 


6 (k*V-'Jj)k* 

mm  rn 


2 ± 


* 

where  oj  = eB/rac  > 0.  On  Che  right  hand  side  of  Eq.  (18),  the  summa- 

c 

tion  over  k includes  also  a summation  over  -k.  That  is  for  each  k in 
the  summation,  there  is  an  equal  term  arising  from  -k.  Multiply  Eq. 
(18)  by  sin  S and  integrate  over  velocity.  The  second  term  on  the 
left  integrates  to  - (u^niy}  The  magnitude  of  the  Integral  of  the  first 
term  is  roughly  p^/L  times  the  integral  of  the  second  term.  Here  is 
the  electron  laraor  radius  and  L in  some  macroscopic  scale  length. 

Thus  as  long  as  p£  « L,  we  can  neglect  the  first  term  on  Che  left  hand 
side  of  Eq.  (18).  The,  using  Eq.  (9)  for  f(V)  on  Che  right  hand  side 
of  Eq.  ( IS) , we  find 
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To  find  the  energy  flux  in  the  reference  frame  moving  with  speed  c 
c Ey/8,  multiply  Eq.  (18)  by  b m(VT  4 V*)V  sin  0 and  Integrate 
over  velocity.  The  same  analysis  as  leading  to  Eq.  (19)  gives 
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The  thermal  energy  flux  corrected  for  fluid  convection  Is  then 
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Thus,  as  long  as  the  plasma  is  unstable  the  energy  flux  is  down  the 
temperature  gradient.  Let  us  note  in  passing  that  if  the  spectrum  is 
one  dimensional  in  the  y direction,  the  arguments  leading  up  to  Eq.  (21) 
shows  there  is  no  instability  generated  heat  flux  in  the  y direction. 

If  one  assumes 


k - 1^/2, 


(22) 


as  is  characteristic  of  ion  acoustic  instabilities,  Eq.  (21)  is  an  ex- 
pression for  the  energy  flux  in  terms  of  the  fluctuating  field  strength. 
Clearly,  it  has  the  same  form  as  the  classical  thermal  energy  flux. 

One  could  in  fact  write  out  an  expression  for  the  anomalous  thermal 
conduction. 
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We  will  close  this  section  by  determining  an  approximate  maximum  value 

for  — for  which  the  magnetic  field  is  Important  in  the  limitation  of 
e 

energy  flux.  An  earlier  study9  has  shown  that  in  an  unmagnetized  plasma, 
an  energy  flux  generates  ion  acoustic  instability  which  inhibits  the 
heat  flux.  There,  it  was  found  that  for  the  unmagnetized  plasma. 
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where  again  we  have  assumed  k ~ 1^/2.  Notice  that  for  the  unmagnetlzed 

I I II 


plasma,  as 


22 
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22 


for 


Increases  shrinks,  where  W Increases  with 
the  magnetized  plasma.  If  is  so  large  that  expressed  in  Eq. 

(1^)  is  smaller  than  W of  Eq.  (21),  then  one  could  reasonably  conclude 


Chat  the  magnetic  field  no  longer  plays  a role  in  limiting  the  heat 
transport.  Instead,  the  energy  transport  is  inhibited  only  by  the  ion 
acoustic  turbulence  and  W is  given  by  Eq.  (2^).  It  is  then  a simple 
calculation  to  show  that  if 


the  magnetic  field  does  not  limit  the  electron  heat  transport. 

The  above  argument  is  even  more  plausible  if  one  makes  the  reason- 
able assumption  that  the  effect  of  the  turbulence  is  to  induce  an 
anomalous  collision  frequency 


In  that  case,  our  expression  for  the  thermal  conduction  coefficient  for 

V2 

IT 

a magnetized  plasma  is  -£■  v jjs-  , while  for  a unmagnetized  plasma,  it 

ce 

is  2 V2/v  . This  is  then  clearly  analogous  with  the  classical  ex- 
e an 

pres sion  for  electron  thermal  conduction12  except  for  numerical  factors 
of  order  unity. 

V APPLICATION  TO  LASER  HYDRO-CODES 

As  shown  in  the  numerical  example  at  the  end  of  Section  III,  mag- 
netized laser  produced  plasmas  will  often  be  unstable.  Since  laser 
fusion  targets  are  generally  designed  by  using  hydrodynamic  codes,22"25 
it  is  essential  to  determine  how  the  presence  of  instability  can  be 
modeled  in  these  codes.  Equation  (23)  gives  an  expression  for  the 
thermal  conduction  in  the  x direction  in  terms  of  fluid  parameters  and 
also  e<c/T#.  It  was  also  stated  in  Section  IV  that  the  instability  has 
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no  effect  on  the  collisionless  heat  flux  In  the  y direction.  To  model 
anomalous  thermal  conduction,  the  first  thing  to  do  is  to  determine  if 
the  plasma  is  unstable  at  a particular  point.  Of  course  ion  Landua 
(and  possibly  collisional)  damping  must  be  figured  into  the  growth  rate. 
If  it  is  stable,  simply  use  classical  thermal  conduction.  If  it  is 
unstable,  the  problem  is  to  determine  ec/T^.  To  do  this,  one  can  make 
contact  with  a great  deal  of  work  on  the  nonlinear  theory  of  ion 
acoustic  instabilities^6-3 4 A general  consensus  seems  to  be  that  some 
sort  of  trapping  or  resonance  broadening  ultimately  limits  the  value  of 
eo/T^  to  somewhere  around  0.1  < eo/T^  <0.2.  Assuming  such  a value  of 
ec/T^,  Eq.  (23)  then  gives  an  approximate  expression  for  the  anomalous 
thermal  conductivity  in  the  unstable  regions.  The  final  thing  to  check 
is  whether  or  not  this  value  of  ec/T^  is  so  large  that  the  turbulence, 
rather  than  the  magnetic  field,  dominates  the  thermal  conduction,  in 
other  words,  if  the  inequality,  Eq.  (25),  is  satisfied.  If  it  is,  then 
instead  of  Eq.  (23),  the  appropriate  thermal  conductivity  is  the  un- 
magnetlzed  value,  either  classical12  or  anomalous.9 

In  summary,  this  paper  has  shown  that  the  thermal  conductivity  in 
laser  produced  plasmas  can  be  anomalously  Increased  due  to  self  genera- 
ted ion  acoustic  turbulence.  Also  it  has  shown  how  this  effect  may  be 
modeled  in  laser  hydro-codes. 
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APPENDIX 


In  order  Co  calculate  whether  a magnetized  plasma  Is  stable,  one's 
first  instinct  is  generally  to  find  a Vlasov  equilibrium  which  describes 
the  unperturbed  situation.  This  is  not  possible  to  do  in  the  plasma 
pictured  in  Fig.  1.  The  problem  is  that  the  pressure  gradient  in  the 
x direction  is  balanced  by  an  electric  field  in  the  x direction  E^,  and 
the  electron  flow  in  the  x direction  is  generated  by  an  electric  field 
in  the  y direction  E = - B.  Thus  neither  P = mV  - eAy/c  nor 

y c y y 

P = mV  - eAx/c,  the  canonical  momenta  in  either  the  y or  x direction, 
x x 

are  constants  of  the  motion.  This  is  quite  different  from  the  case  of 
a magnetically  confined  plasma,  where  there  is  no  flow  velocity  in  the 
x direction  so  that  P^  is  a constant  of  the  motion.  This  allows  dis- 
tributions with  x variations  in  density  and  temperature  to  be 
constructed  by  exploiting  the  constancy  of  P^. 

In  the  laser  produced  plasma  shown  in  Fig.  1 one  can  see  in  another 
ray  that  there  is  no  relevant  Vlasov  equilibrium.  To  do  so,  notice 
that  a Vlasov  equilibrium  implies  no  dissipation.  Then,  if  the  dis- 
tribution function  is  nearly  Maxwellian,  ideal  fluid  equations  should 
apply.  That  is,  as  the  fluid  flows  into  lower  density  regions,  it 
should  cool  adiabaticaliy , as  it  expands  or  T ~ n2/3.  The  density  and 
temperature  gradients  then  could  not  possibly  be  in  opposite  directions 
as  shown  in  Fig.  1. 

The  reason  the  density  and  temperature  gradients  are  opposite  for 
the  laser  produced  plasma  is  that  the  fluid  is  not  ideal;  its  equilib- 
rium properties  are  dominated  by  electron  thermal  conductivity. 
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a b 

Fig.  2 - (a)  The  lines  of  constant  f(VxVy)  in  the  frame  moving 
with  speed  cEy/B.  Dotted  lines  are  for  the  absence  of  insta- 
bility, solid  lines  are  for  the  presence  of  instability,  (b) 
Lines  of  constant  f(VxVy)  in  the  frame  with  u = 0 with  the 
presence  of  instability. 
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